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Abstract 

We investigate saturation effects in susceptible-infected-susceptible (SIS) models of the spread 
of epidemics in heterogeneous populations. The structure of interactions in the population is 
represented by networks with connectivity distribution P(k), including scale- free (SF) networks 
with power law distributions P{k) ~ k~~' . Considering cases where the transmission of infection 
between nodes depends on their connectivity, we introduce a saturation function C(k) which reduces 
the infection transmission rate A across an edge going from a node with high connectivity k. A mean 
field approximation with the neglect of degree-degree correlation then leads to a finite threshold 
A c > for SF networks with 2 < 7 < 3. We also find, in this approximation, the fraction of infected 
individuals among those with degree k for A close to A c . We investigate via computer simulation 
the contact process on a heterogeneous regular lattice and compare the results with those obtained 
from mean field theory with and without neglect of degree-degree correlations. 

PACS numbers: 89.75.-k,87.23.Ge,05.70.Ln 
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I. INTRODUCTION 



Many social, biological, and physical systems can be modeled as networks, i.e., connected 
graphs with at most a single edge between nodes: nodes represent entities and edges represent 
interaction pathways among the nodes The connectivity pattern in these networks 

encode information about the structure of the system 0, U, |||. An important and much 
studied feature of these networks is their degree distribution P{k), where P(k) is the fraction 
of nodes of the network that have k connections to other nodes. It was found that many 
interesting networks such as the internet p and the patterns of human sexual contacts 
Ml are very heterogeneous with approximately " scale- free" (SF) degree distribution: i.e., 

nnnu 

P{k) ~ k 7 (power law distribution) with 2 < 7 < 3 |a 12|, |8|, |9| . The study of epidemics in 
heterogeneous networks is therefore of practical importance for the control of the spread of 
cyber viruses and biological epidemics. 

The mathematical studies of epidemics on the other hand often make the assumption of 
a homogeneous population Q, Q] . This means that any infective individual is an equally 
likely source for the further transmission of the disease to other members of the population 
with whom that individual is in contact and vice versa. The simplest epidemiological model 



of that kind is the SIS model 
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111 ]. In the SIS model, individuals can only exist in two 



discrete states; healthy (but susceptible) or infected. The disease processes are specified as 
follows: Infected individuals become susceptible (healthy) at rate 5, independently of their 
environment. We shall choose time units in which 5 = 1. Susceptible individuals become 
infected at a rate A multiplied by the number of infected neighbors, i.e., infected nodes to 
which they are connected by an edge. When the web of interactions between individuals 
is taken to be a reg ular lattice the stochastic process describing this system is the Harris' 



contact process 



121. 



Epidemic behavior in these homogeneous networks, where each node has z neighbors, 
P{k) = 5k,zi show a "phase transition" as the rate A at which an infected individual infects 
a susceptible neighbor, is changed; i.e., there exists a critical value A = A c > below which 
the only stationary state is a disease-free state (or absorbing phase) and above which there is 
an endemic infected state (or active phase). This can be proven rigorously for the stochastic 
contact process on a regular lattice, an infinite homogeneous network and is inherited by 
mean field models based on this process. The mean field critical value, X^ F , is proportional 
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to z -1 , the inverse of the number of "interacting neighbors" Q. I"!^. This mean field \^ IF 
is smaller than that for the contact process. The latter depends not only on the number of 
neighbors but also on the topology of the lattice (see later) 12. Il3|. 

An interesting question then is how to extend these models, which correspond to net- 
works with homogeneous connectivity, to real world situations where the number of contacts 
varies greatly from one node to another [3, In such heterogeneous net- 

works, each node has a statistically significant probability of having a very large number of 
connections compared to the aver age connectivity of the network. The mean field version 

of this problem was studied in 0,0 0,0, HQ where it was shown that the epidemic 

threshold decreases with increasing second moment of the connectivity distribution. As a 
result epidemic processes in infinite SF networks with diverging second moment, 7 < 3, are 
believed not to possess any epidemic threshold below which the infection cannot produce an 
epidemic outbreak or an endemic infected state 

HI 

The absence of an epidemic threshold in SF networks makes them very vulnerable. This 
remains true even if one takes into account the finite size of real systems which of course 
always have finite second moments. In general the epidemic threshold for a heterogeneous 
network is much smaller than for a homogeneous network with the same average number of 
contacts |l9j. The presence of " assortative" or " dissortative" two-point degree correlation in 



SF networks with 2 < 7 < 3 does not appear to alter the absence of epidemic threshold 
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In these analyses all edges are treated in the same way. There are however many situations 
where there are differences between the "strength" of different edges. We investigate here 
cases where the assigned weight of an edge between two nodes depends on the connectivity 
of these nodes. We consider in particular saturation effects due, for example, to the fact that 
disease transmission requires a certain amount of contact time or "finite time commitment" 
from both individuals in contact. This has the effect of lowering the effective connectivity 
for highly connected nodes and thereby decreasing the importance of the "heavy tails" in 
P{k). Such saturation effects then lead to a finite threshold even when the second moment 
of P(k) diverges. An example in which saturation effects, due to temporal limitation of 
interactions, play an important role is Holling's "the principle of time budget" in behavioral 
ecology |23- 

Using mean field approximations appropriate for heterogeneous network models ll^ 
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of epidemics we calculate the critical value A c for different saturation patterns. We also find 
in this mean field approximation, which neglects degree-degree correlations between different 
nodes, the behavior of the endemic prevalence pk, i.e., the fraction of infected nodes of degree 
k, for A close to A c . They all have the same behavior for A ~ A c , pk ~ Ak(X — X c ) , with 
increasing with k; (3 = 1 when the third moment is finite. This dependence on k is missed 
by the homogeneous approximation of the contact network. 

We then investigate, via numerical simulations, the behavior of the stochastic contact 
process on a regular lattice consisting of nodes with two different degrees. The results 
are compared with mean field approximations with and without neglect of degree-degree 
correlations. 



II. MEAN FIELD SIS MODEL WITH SATURATION 

The mean field theory for the contact process, obtained by neglecting correlations between 
different nodes, is described by an equation for the density of infected nodes pit), present 
at time t, which can be written as [13] 

^ = -p(t) + \zp(t)(l-p(t)) (1) 

where z is the coordination number. Solving Eq. yields p(t) = (\z — l)p(0)e < - Az_1 ^/[A-2 — 
1 + Azp(0)(e( Az-1 )* — 1)]. The steady state solution of Eq.(l), obtained as t — ► oo, has an 
epidemic threshold X^ 1F = 1/z. For A > A* /i? , any initial infection spreads and becomes 
persistent with stationary total prevalence level p = (A — \^ IF )/\. Below the threshold(A < 
), the initial infection dies out exponentially fast. A similar transition occurs for the 
stochastic contact process on a regular lattice with edges between nearest neighbors. The 
critical values A c for the lattice Z d are A c ~ 1.6489 in one-dimension (d = 1, z = 2), 
A c ~ 0.4122 in two- dimension (d = 2, z — 4), etc. z\ c approaches z\^ IF = 1 as d — ► oo [13]. 

Consider now a general network with degree distribution P(k). Let pk{t) be the fraction 
of the nodes with degree k which are infected at time t. We define \C(k, I) as the effective 
transmission, or infectivity, rate across an edge going from a node with degree A; to a node 
with degree I: C(k,l) = 1 in the absence of saturation effects. 

nn 

The mean field equation of the contact process on this network |15l . Il6| which ignores 



correlations between the states of the nodes, yields the following set of differential equations 
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for p k (t) 

^- = - Pk (t) + XQ k (t)(l-p k (t)) (2) 

We can interpret XQk{t) as the effective transmission rate of infection to an uninfected node 
of degree k by all infected nodes with which it is in contact via any of its k edges, 

Q k (t) = kJ2P(l\k)C(k,l) Pl (t) (3) 
i 

Here P(l\k) is the probability that an edge emerging from a node of degree k has its other 
end at a node with degree I and C(k,l) is the effective strength of such a bond. We shall 
now assume further that 

P(l\k) = (4) 

with z = Ylk kP(k), i.e. random attachment (no degree-degree correlation). 

To specify C(k, I) in Eq. (j3*|) . we make the simplifying assumption that an individual with 
k contacts spends equal time with each neighbor. The effective strength of an edge is then 
given by a product of ratios of effective connectivity to total connectivity of each node in 
contact, 

C(k, I) = C{k) ° {l) (5) 

with C(k,l) = 1 corresponding to uniform bond strength. 
Eq. (j2J can now be written in the form 

^ = - Pk (t) + A(l - Pk (t))C(k)e({p(t)}) (6) 

where 

emm = E^(*)C(*)/»M (7) 

Multiplying Eq. (jHJ) by c ( fc ) p ( fc ) anc l summing over k yields 

^ = -e(t) + ^ E P(k)c\k)(i - Pk (t)) (8) 

k 

where Q(t) is shorthand for 0({p(t)}). Given P{k) and C(k), Eqs. ©, © and © form a 
closed set of nonlinear differential equations for the Pk(t) which can be solved in principle for 
any given initial values {pj(0)}. They reduce to a single equation, Eq. ((TJ, when P(k) = 5k, z 
and p z = p. For a general P(k) the number of variables and equations are infinite. 
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We are interested primarily in finding stationary solutions pk, < pk < 1, in which not 
all pk = 0. (There is of course always one solution corresponding to the uninfected state 
Pk(t) = Pfc(0) = for all k.) We write the stationary version of Eq. (0) in the form 

C{k)XQ 

Pk ~ i + c(k)\e (9) 

Multiplying Eq. (jHJ) by p ^ c ^ an d summing over k we get 



where f(x) is a monotone decreasing function of x. Eq. (fTU|) will have a (unique) solution 
0(A) different from zero if and only if A > A c , the epidemic threshold, 

1 z 

Xc = W) = WW) 

For diverging second moment, (C 2 (k)) = oo, A c = 0. For a regular lattice with P(k) = Sk >z 
and C{k) = k, we recover the usual mean field result, A c = -. 

Once we have found 0(A) > from Eq. (fTU|) we then get pfc(A) > directly from Eq. (JHl) 
for all k for which C(k) > 0. To get the behavior of the /5fc(A) for A j, A c we expand the right 
side of Eq. flTUj) in small AO. This yields 

T = r - V S P ^)C\k) + 0((AO) 2 ) (12) 

c k 

Solving for O with finite (C 3 (k)), we get 

and we obtain from Eq. (JHJ), 

p k ~ A C C(A;)A(A - A c ) (14) 



p = ^ P(A;)p fc ~ A C (C(A;))A(A - A c ) (15) 

k 

where A = Af#M ' 

When (C 3 (k)) = oo both O and p are not differentiable as A approaches A c from above. 
For a SF network with C{k) = k, (C 3 (k)) is finite only for 7 > 4 where p ~ (A — A c ). 
However for SF network with connectivity saturation, the range of 7 where (C 3 (k)) is finite 
can include all cases with 7 > 2 so that z is finite. The case where (C 3 (k)) is infinite is 
discussed in the Appendix. 



III. EPIDEMIC THRESHOLD IN SF NETWORK WITH CONNECTIVITY SAT- 
URATION 



To see how connectivity saturation modifies the behavior of epidemics, we consider two 
different types of saturating functions C(k), 

Cj(k) = < k iik< kmax (16) 

k if k > k 

""max " "> C "'max 

"'max \ ™ & 

where k is total connectivity of a node and k max is a parameter. 

We replace the sum in Eq. fllOl) by an integral over /c and carry out calculations of A c for SF 
networks with P(k) = (7 — l)m 7_1 A;~ 7 for k > m, 7 > 2. After elementary integration we can 
obtain the epidemic thresholds A c by using Eq. (fTTj) and the second order moment (C 2 (k)). 
They are plotted against 7 in Fig. ^ This figure presents a phase diagram consisting of two 
phases: a disease-free state below each epidemic threshold curve and an endemic infection 
state above each curve. Note that A c diverges as 7 — > 2 as can be seen from divergence of z 
in Eq.(ll) when (C 2 (k)) is finite. This can be understood by noting that as the number of 
edges increases the effective infection rate for any node decreases. Note also that A c = in 
the absence of saturation for 2 < 7 < 3 because of the divergence of (C 2 (k)). 

In the inset of Fig. ^ we compare the epidemic thresholds A c of SF networks with satura- 
tion with A^ omo , obtained from a homogeneous network with coordination number z — (k). 

The dependence of the epidemic threshold on k max is plotted in Fig. El When k max 
is finite, the epidemic threshold A c (7) is non-zero for 7 > 2 and as k max increases A c (7) 
decreases. When k max = 00, the epidemic threshold of SF network without saturation is 
recovered. 

The stationary total prevalence p for the SF network with and without saturation are 
plotted in Fig. El The stationary total prevalence p(A) with saturation is smaller than that 
without saturation for all A > 0. This is because saturation reduces the effective transmission 
rate of infection to an uninfected node across an edge going from an infected node with high 
connectivity (see Eq. Q and (JSJ)). 



7 



IV. THE CONTACT PROCESS AND MEAN-FIELD APPROXIMATION IN A 
HETEROGENEOUS REGULAR LATTICE 



To investigate the effect of heterogeneity and degree- degree correlations we investigated 
the SIS model on the "face-centered" square (FCS) lattice with two types of nodes. Type 
A nodes, which connect to both nearest and next nearest neighbor sites, have connectivity 
kA = 8 and type B, which connect only to nearest neighbor sites, have ks = 4. For this 
system P(8|4) = 1 and P(4|8) = P(8|8) = 1/2, and P{k y \k a ) = otherwise: see Fig.H 

We carried out computer simulations on this model with no saturation. The critical point 
A c and critical exponents, 5, v» and vj_, were obtained by using the dynamical Monte Carlo 
method j^. As expected there is a critical point A c = 0.23(6). This value of A c is closer to 
A c = 0.18(1) for the homogeneous regular lattice with z = 8. than to A c = 0.412 jl3j] for the 
square lattice with z — 4. The homogeneous regular lattice with z = 8 is the square lattice 
with both nearest neighbor (NN) and next NN bonds (see Fig. The critical exponents 
appear to be the same as for the square lattice as expected from universality considerations. 
The phase diagram is plotted in Fig. El 

We also computed Pa and ps for A > A c . The results are plotted in Fig. |BJ We can see 
that nodes with higher connectivity are more infected than those with less connectivity at 
A close to A c . 

The mean field equations of the SIS on this lattice with the exact P{k 1 \k a ) are, 

dp A 



dt 
dp B 
dt 

The steady state solutions are given 



~Pa + 4A(1 — Pa){pa + Pb) (18) 
-p B + A\p A {l- p B ) (19) 



PA = ^P^y PB = 1 -^7CT (20) 

with A c = (-1 + V5)/8. We call this "MF1". 

The mean field equations with "no degree-degree correlation" approximation, "MF2" , can 
be obtained in a similar manner to that given in Eq.(6)-(7). This corresponds to putting 
P(4|4) = P(4|8) = 1/3 and P(8|4) = P(8|8) = 2/3. 

The results from simulation and the two mean field approximations are compared in 
Fig. and El and in Table. HI Close to the critical point, we can approximate the ratio, 
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Pb/pa, by neglecting nonlinear term in Eq. (19): Pb/pa ~ 4A C P(8|4)/(1-4A C P(4|4)). In this 
heterogeneous lattice with only two types of nodes, "MF2" with no degree-degree correlation 
gives the same A c but not as good values for p A and p B as "MF1" with the exact degree- 
degree correlation. 

V. CONCLUDING REMARKS 

In this manuscript we considered the SIS epidemic model on heterogeneous networks with 
saturation. This made the epidemic thresholds finite for SF networks with 2 < 7 < 3. We 
also investigated via computer simulation the stochastic contact process on a heterogeneous 
regular lattice and compared the results with those obtained from mean field theory with and 
without neglect of degree-degree correlations. Our considerations extend naturally to other 
types of heterogeneous networks in which the effective strength of an edge depends on the 
degrees of the nodes which it connects. Thus in considering the spread of computer viruses 
on the internet, effects similar to saturation might arise from nodes with high connectivity 
having higher "firewalls" around them. 
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APPENDIX A: CRITICAL BEHAVIOR OF p WHEN (C 3 (k)} = 00 

The critical behavior of the steady state 6 and p in the presence of saturation can be 
evaluated by using the third order moment (C 3 (k)). When k max < 00, both G and p can 
be expanded in a power series close to the critical point because of z, (Cf(k)) and (Cfj(k)) 
being all finite for 7 > 2. 

When k m ax — *■ 00, (C 3 (k)) may diverge and as a result the expression of Eq.(13)-(15) are 
not valid any more. Let us introduce the limiting case of Cn{k) when k max = 00, 

C ni (k)=k? for^<p<l (Al) 
9 



Then, {Cjjj(k)) = for 7 > np + 1 while (Cf n (k)} = 00 for 2 < 7 < np + 1. 

When 7 > 3p + 1, i.e., {Cj II (k)} < 00, 6 and p at A close to A c are equivalent to those in 
Eq.(13)-(15). 

However when 2 < 7 < 3p + 1, i.e., (C] n (k)) = 00, p and are treated in different 
way and their behaviors close to the critical point are presented as follows: For Cjn(k) the 
steady state Q Cl11 and p Cni can be given, 



Q Cin = ( 7 ^^ ^[M-M.-J-] (A2) 
2(7 - p - 1) A0m p 



for non- integer values of 77 = (7 — l)/p and for 7 > 2. After Taylor-expansion in the limit 
of small we obtain (i) p Cl " ~ Xto-V/P-v) for 2 < 7 < 2p + 1 (ii) p^" ~ A x /(^ 2 ) for 
2p + 1 < 7 < 3p + 1 
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TABLE I: The critical point A c , the ratios ipA(B) an d 4>, ipA(B) = H m X—f\ c PA \ B \ C an d 4> = 
lim\^\ c ^, of the SIS epidemics on the face-centered square lattice. "Sim" refers to simula- 
tion data without saturation(i.e., C(k) = k) while "MF1" and "MF2" refers to mean field analysis 
with exact P{k 1 \k a ) and approximate one using Eq.(4), respectively. 

Sim 0.23(6) oo oo 0.7 

MF1 0.15 7.24 4.5 0.62 

MF2 0.15 7.4 3.7 0.5 




FIG. 1: Epidemic thresholds A c as a function of 7. Thin and thick solid lines are drawn for 
epidemic thresholds from SF networks with saturation type Ci and C/j, respectively. Fat solid 
line represents the epidemic threshold from SF network without saturation. The dashed line is the 
epidemic threshold for a homogeneous network with coordination number z = (k). Inset: Ratios 
of the epidemic thresholds for the SF network with saturation and without saturation to that of 
the homogeneous network. k max = 100, m = 5 and p = 1 are used. 
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FIG. 2: Dependence of epidemic threshold from SF network with saturation function Cu on the 
cut-off connectivity k max . From top to bottom k max = 10, 10 2 , 10 3 , 10 4 and oo. Here p = 1 for 
C„(k). 



13 




FIG. 3: The steady state p(A) and 6(A) of the SIS epidemics on the SF network for 7 = 3 as a 
function of infectivity A. Thick solid(thick dashed) lines indicates the steady state p Cl (6 C/ ) with 
Type I saturation. Thin lines are for those without saturation. The critical point in this particular 
case is given: X^ 1 = 0.0573 for Type I saturation and A c = for no saturation. Here k max = 100 
and m = 5 are used. 
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(b) 

FIG. 4: The Topology of lattices, (a) The face-centered square (FCS) lattice. Nodes with type A 
have 8 degrees while nodes with type B have 4 degrees, (b) A homogeneous network with z = 8. 
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FIG. 5: The steady state prevalence p of SIS in the FCS lattice without saturation as a function 
of infectivity rate A. Pa(pb) represents the fraction of infected nodes of type A and type B, 
respectively and p = {pA + /5b)/2. Simulation ("sim") data are shown with circles while mean 
field results, where "MF1" with exact P{k 1 \k a ) and "MF2" with approximate one using Eq.(4), 
are given with solid and dashed lines, respectively. Inset: Simulation results of the contact process 
in the FCS lattice (circles), in the square lattice (triangles) with z=4 and in the square lattice 
with both nearest and next nearest neighbor interactions (squares) with z=8. All simulations were 
performed on lattices of size L 2 = 100 2 with random initial distribution of infected nodes. 
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FIG. 6: Ratio of pB to pA in the FCS lattice as a function of infectivity rate A. The notation 
and symbols are the same as in Fig.4. The arrows point to the critical A c corresponding to each 
scheme. 
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